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Abstract 

Blunt-forebody pressure data are used to study the 
behavior of the NASA Dry den Flight Research Center 
flush airdata sensing (FADS) pressure model and 
solution algorithm. The model relates surface pressure 
measurements to the airdata state. Spliced from the 
potential flow solution for uniform flow over a sphere 
and the modified Newtonian impact theory, the model is 
shown to apply to a wide range of blunt-forebody shapes 
and Mach numbers. Calibrations of a sphere, spherical 
cones, a Rankine half-body, and the F-14, F/A-18, X-33, 
X-34, and X-38 configurations are shown. The three 
calibration parameters are well-behaved from Mach 
0.25 to Mach 5.0, an angle-of-attack range extending to 
greater than 30°, and an angle-of-sideslip range 
extending to greater than 15°. Contrary to the sharp 
calibration changes found on traditional pitot-static 
systems at transonic speeds, the FADS calibrations are 
smooth, monotonic functions of Mach number and 
effective angles of attack and sideslip. Because the 
FADS calibration is sensitive to pressure port location, 
detailed measurements of the actual pressure port 
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locations on the flight vehicle are required and the wind- 
tunnel calibration model should have pressure ports in 
similar locations. The procedure for calibrating a FADS 
system is outlined. 
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Nomenclature 

angle-of-sideslip solution parameters 
angle-of-attack solution parameters 
angle-of-sideslip solution parameters 
pressure coefficient, (p - p^)/q x 
impact pressure coefficient, (p - p oa )/q c 
diameter 

electronic scanning pressure 

., _ . major axis 

ellipsoid fineness ratio, —d -— 

minor axis 

flush airdata sensing 

FADS pressure model 

High Alpha Research Vehicle 

Lockheed Martin Vought Systems 
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National Technical Systems 
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total pressure, lbf/ft - 
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Subscripts 

e 

i,j,k 

max 

ref 


total pressure behind a normal shock 
wave, lbf/ft 2 

individual port weightings 

impact pressure, lbf/ft 2 
(subsonic q c = p t - p m ; 
supersonic q c = p - p m ) 

free-stream dynamic pressure, lbf/ft 2 

weighting matrix 

Reynolds number 

Systems Research Aircraft 

Unitary Plan Wind Tunnel 

velocity 

Taylor series parameter 

Cartesian coordinates 

FADS pressure model parameters 

angle of attack, deg 

effective, or local, angle of attack, deg 

angle of sideslip, deg 

effective, or local, angle of sideslip, deg 

shape and compressibility parameter 

e component caused by Mach number 

maximum acceptable Mach error 

angle-of-attack flow correction angle, deg 

angle-of-sideslip flow correction angle, 
deg 

ratio of specific heats (1.4 for air) 

pressure difference, = p i - p -, and so 
forth 

local flow incidence angle, deg 
cone angle of a pressure port, deg 
clock angle of a pressure port, deg 
standard deviation 


effective (or local) angle 
port indices 
maximum value 
reference value 


true 

true value 

OO 

free-stream value 

Superscripts 

estimate 

(m) 

iteration number 

T 

matrix transpose 

-1 

matrix inverse 


Introduction 

The accurate measurement of airdata is critical to the 
flight control, guidance, and postflight analysis of most 
modem atmospheric flight vehicles. The entire airdata 
state can be described by five parameters: Mach number 
(M), angle of attack (a), angle of sideslip (P), either 
pressure altitude or free-stream static pressure (p^), 
and the true airspeed. Using these five parameters, all 
other airdata parameters of interest can be calculated. 

Historically, airdata have been measured through the 
use of intrusive booms or probes that penetrate the flow 
away from the influence of the vehicle body to measure 
total and static pressure, angle of attack, angle of 
sideslip, and free-stream temperature. However, 
specialized requirements of advanced vehicles make 
using intrusive conventional airdata measurement 
systems 1 highly undesirable. Advanced vehicles include 
hypersonic vehicles, 2 for which the heating environment 
is extremely hostile; stealth vehicles, which require a 
minimal radar cross section; and high-performance 
vehicles 2 that operate beyond the poststall flight regime. 
For example, the presence of an intrusive airdata 
noseboom on the X-31 aircraft induced unsteadiness in 
the primary forebody vortices, causing unwanted lateral 
bistabilities in the poststall flight regime. 4,5 

The flush airdata sensing (FADS) system concept, in 
which abdata are inferred from nonintrusive surface 
pressure measurements, was developed to circumvent 
many of the difficulties with intrusive airdata systems. 
This method does not require probing of the local flow 
field to compute abdata parameters, and allows for 
airdata measurements in flight applications where 
probe-based systems are no longer viable. 

The FADS concept relies on a mathematical model 
that relates measured forebody surface pressures to the 
airdata state. The current method” uses a nonlinear 
pressure model and a nonlinear solution algorithm to 
solve for the airdata state from the measured pressures. 
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Although the FADS pressure model (FPM) has been 
validated for a high-performance fighter aircraft 

a 

forebody and several simple shapes (for example, 
cylinder and sphere), the suitability of the model for a 
wide variety of vehicle shapes has not been established. 

The primary objective of this study is to gather a set 
of forebody pressure data for a series of blunt-body 
shapes to generally characterize the behavior of the 
FPM calibration. A set of guidelines is developed to aid 
the design of future blunt-body FADS systems. 

The forebody pressure data have been obtained from 
five sources: a wind-tunnel test of three simple blunt- 
forebody configurations; a collection of previously 
published wind-tunnel data for generic blunt-body 
shapes; wind-tunnel tests of the X-33, X-34, and X-38 
hypersonic reentry vehicles; exact potential flow 
solutions for uniform flow over a cylinder, sphere, and 
an arbitrarily-oriented ellipsoid; and flight data from the 
F-14 and F/A-18 aircraft. The nature of the calibration 
for each of the configurations is evaluated as a function 
of Mach number, angle of attack, and angle of sideslip. 
Behavior of the mathematical model is also evaluated at 
Mach numbers much greater than have been previously 
tested (a maximum of Mach 5.0), which will directly 
benefit hypersonic FADS applications that are currently 
in development (the X-33, X-34, and X-38 vehicles). In 
addition, the adaptability of the FPM to nonspherical 
forebody shapes is tested. The results presented 
demonstrate that the current FADS methodology is 
universally applicable to blunt-nosed configurations. 

A patent has been filed on this NASA invention. 
Equations (10)—(16) are included in this patent. Note 
that use of trade names or names of manufacturers in 
this document does not constitute an official 
endorsement of such products or manufacturers, either 
expressed or implied, by the National Aeronautics and 
Space Administration. 

Background and Historical Perspective 

The original flush airdata system prototype was 
developed for the X-15 program and used a steerable- 
ball nose that attempted to nullify upper and lower 
pressure measurements to determine the stagnation 
point and the local angle of attack. 6 The mechanical 
design of this system was cumbersome, and the analog¬ 
steering concept was abandoned after the X-15 program 
ended. 

A subsequent approach, the Space Shuttle entry 

n 

airdata system, was developed at the NASA Langley 


Research Center (Hampton, Virginia) for the Space 
Shuttle program. This system used a matrix of fixed 
static-pressure orifices on the orbiter nosecap. The 
technique later was adapted to aeronautical applications, 
and several demonstration programs were performed in 
the early 1980’s at the NASA Dryden Flight Research 

O Q 

Center (Edwards, California). For these early 
programs, the sensed pressures were related to the 
desired airdata parameters using a nonphysical 
mapping. These tests verified the feasibility of the fixed- 
orifice concept but did not demonstrate real time- 
capable algorithms for estimating the airdata from the 
pressure measurements. 

The first estimation algorithms capable of real-time 
operation were developed at NASA Dryden during the 
late 1980’s for the F/A-18 High Alpha Research Vehicle 
(HARV) program. The HARV flight tests also 
demonstrated that the measurement range of FADS 
systems could be extended to angles of attack greater 
than 60°. The computations were performed postflight 
using pressure data telemetered to the ground. The 
estimation algorithms developed for the HARV program 
recently were demonstrated in a real-time flight 
environment on the NASA Dryden F/A-18 Systems 
Research Aircraft (SRA). 10 

For the HARV and SRA FADS systems, surface 
pressure measurements were related to the desired 
airdata parameters using a calibrated aerodynamic 
model. The model was derived from a combination of 
potential flow over a sphere” and modified Newtonian 
flow theory. 12 This model, which was calibrated with 
flight data, captured the salient features of the local flow 
but was simple enough to be invertable in real time. 
Nonlinear regression was used to invert the 
aerodynamic model. In this algorithm, all surface 
pressure measurements were simultaneously used to 
infer the airdata by linearizing the equations around the 
result of the previous data frame. 

The nonlinear regression algorithm developed for the 
F/A-18 flight test programs exhibited problems with 
algorithm stability in the transonic and supersonic flight 
regimes and in the presence of undetected sensor 
failures. These stability problems required ad hoc 
software patches to artificially aid the stability by 
deleting failed or degraded port measurements from the 
estimation algorithm. These ad hoc additions to the code 
were computationally cumbersome and did not 
universally stabilize the algorithm for all flight 
regimes. 10 The instability problems exhibited by the 
nonlinear regression algorithm precluded its use in a 
real-time feedback mode. 
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To avoid problems encountered using the nonlinear 
regression algorithm, a new solution algorithm 2 has 
been developed for the X-33, X-34, and X-38 
demonstration vehicles. These transatmospheric 
vehicles generally require airdata measurements for 
flight-critical tasks (such as inertial guidance 
stabilization and inner- and outer-loop flight control) 
and for terminal-area energy management. This solution 
algorithm, superior to the nonlinear regression 
algorithm, has been developed by taking strategic 
combinations of three sensor readings to analytically 
decouple the angle-of-attack and angle-of-sideslip 
computations from the Mach number and static pressure 
calculations. 2 This innovation allows development of an 
estimation algorithm whose solution speed is superior to 
the nonlinear regression algorithm and whose solution 
stability can be analytically assured for a given port 
arrangement. 

Flush Airdata Sensing Pressure Model 

As previously mentioned, the FPM is used to relate 
the observed surface pressures to the desired airdata 
state. The model is postulated as a compromise between 
a simple potential flow model on a sphere, 11 


c p ( 0) = 


p(V)-p. 


1 9 ■ 2 o 

= 1 - - sin 9 

4 


( 1 ) 


5 9 2„ 

— + -cos 9, 
4 4 


and modified Newtonian flow theory 12 for blunt objects 
in hypersonic flow, 


Pt 2 ~P - 9 

Cf9) = —-cos" 9 = C cos" 9 , (2) 

P a Pm ax 

1 oo 

where C p is the pressure coefficient at a point on the 
forebody surface and 9 is the incidence angle between 
the surface normal at the pressure port and the local 
velocity vector. The incidence angle is geometrically 
related to the effective (local) angle of attack, a g , and 
effective (local) angle of sideslip, P g , by 10 

cos(9 ( ) = cos(a g )cos(P g )cos(X ( ) 

+ sin(P c )sin(<|> ( )sin(?i ( ) (3) 

+ sin(a e )cos(P g )cos((|) ( )sin(?i.), 

where i is a port index, hr equation (3), cp and A. are the 
local surface coordinates referred to as the clock and 
cone angles, respectively. The clock angle is measured 


clockwise around the longitudinal axis looking aft 
starting from the bottom. Cone angle is the total angle 
the normal to the surface makes with respect to the 
longitudinal axis of the forebody. Figure 1 shows the 
clock and cone angle definitions. 




Figure 1. Clock and cone angle definitions. 


The effective angles of attack and sideslip, r/. f , and 
P g , are the flow direction angles as locally sensed on the 
blunt forebody. These angles have been deflected from 
the free-stream flow direction by expansion around the 
forebody and by upwash and sidewash induced by the 
rest of the vehicle. Tire true free-stream flow incidence 
angles are related to the effective angles by a set of 
calibrations" of the forms 

§ a [M 00 ,a g ] = a e -a true (4) 

and 

SPIM^P,] = P e -P, r „ e . (5) 

Tire calibrations must be empirically determined using a 
reference airdata set obtained either in a wind tunnel or 
from flight data. Fortunately, the FADS wind-tunnel 
calibration tests can be easily accomplished in parallel 
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with traditional force and moment testing. Methods for 
determining these calibrations from a reference set are 
described in detail in the “Flush Airdata Sensing 
Pressure Model Calibration Procedure" section. 

Equations (1) and (2) can be spliced together to give a 
flow model that is applicable to a large Mach number 
range by the parameterization 

C p (0) = X(M oo ,a g ,P g ) + Y(M oo ,a e .P g )cos 2 0 , (6) 

where the coefficients X and Y are a function of Mach 
number and the effective flow angles. At 0° incidence 
angle, equation (6) must satisfy subsonic flow 
conditions: 


C„(0) = 


Pt~P oo q c 


7oo <7oo 

or supersonic flow conditions: 

Pu-Po o q c 

C_(0) = -2- = — . 

P <7oo 7oo 


(7(a)) 


(7(b)) 


By setting 


X(M TO ,a e , P g ) = —e(M oo ,a g ,P e ) (8) 

i oo 


and 


Y(M oo ,a e ,P g ) = ^[l-e(M co ,a g ,p g )], (9) 

the constraints of equation (7) are satisfied for all Mach 
numbers and the general FPM results: 


effective flow vector and the surface normal, and a 
calibration parameter, e. Because flow incidence angles 
a e and p r are imbedded in Band e, equation (10) is 
inherently nonlinear. Usually five or more pressure 
measurements are used to solve for the minimum airdata 
state ( q c , ;) M , a, p) , 2, 10 although the airdata can be 
calculated with as few as four. The minimum airdata 
state is sufficient to calculate all other important airdata 
parameters (pressure altitude, Mach number, equivalent 
airspeed, and so forth) except true airspeed, which 
requires a measurement or estimate of the free-stream 
temperature in order to calculate the speed of sound. 

As was mentioned in the Introduction, the nonlinear 
regression algorithm that was applied in references 3 
and 10 had problems with algorithm stability. A 
preferable method for estimating the airdata parameters 
is described in the following subsections. A detailed 
derivation and stability analysis has been published,* 

The “Triples" Algorithm 

By taking strategic combinations of three surface 
sensor readings (“triples”), q c , p M , and e are eliminated 
from equation (10). The resulting equation is 

Pi~ Pi 
Pj-Pk 

2 2 2 2 

q c { cos 0 j + £sm 0,} + p„ - [q c { cos 0 ( + Esin 0.} + p„J 

q c { cos 2 e ( . + esin 2 e ; .} + p„ - [q c { cos 2 0 t + Esin 2 0*} + p J 

... , 2 „ , ... , 2 „ . (11(a)) 

{(1 -e)cos 0 ; + e}-{(l -e)cos 0 ( + e} 

{(1 - e)cos 2 0 ( -+ e } - {(1 - e)cos 2 0 <: + e} 

2 „ 2 „ 

cos - cos 8 ■ 

_ 2 2 ' 

cos 0• - cos 0£ 

Rearranging equation (11 (a)) yields 


p(Q) = 

q c [cos 0(oc e ,P g ) + e(M oo ,a g ,P g )sm 2 0(a g ,P g )] (10) 

+ p . 

I oo 

In equations (7)-(10), q^ is the free-stream dynamic 
pressure and q c is the impact pressure. 

Flush Airdata Sensing Solution Algorithm 

Equation (10) states that the local surface pressure at 
any point on a spherical (blunt) forebody is a function of 
the free-stream flight conditions, the angle between the 


F* cos 2 + F ( . cos 2 0^ + r kj cos 2 0 ( = 0 (11(b)) 

where 

T ik = ( Pi~Pk ); 

T j> = ( Pj~ p > (ihc)) 

r kj = ( Pk-Pj ); 

and pj, pj, and p k are the pressures used in the data 
“triple.” Although equation (11(a)) is still nonlinear, e 
has been removed, thus decoupling the local flow angles 
from the calibration parameter. 
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Angle-of-Attack Estimator 

The effective angle of attack, a e , can be further 
decoupled from P f by using only pressures aligned 
along a vertical meridian (where <|) = 0° or 180°). In 
this geometry arrangement, terms related to (t r are 
eliminated from equation (11(b)). The result is a 
quadratic expression in tan(a ) of the form 


T ik sin2 X j + F j, X k + F kj sin2 \ 


. 2 
tan a 


A ' = ^■*b/ + r. i b t 2 + r t .b l 2 


B = 1 a j h j + O' a A b A- + T kj a i h i 


C = ^■*a J 2 + r y ,a/ + r ty .a 1 . 2 


and 


+ 2 


T ik cos^sinXyCOSt^ + IF cos X k sin X k cos(| 


( 12 ) 


a {ijk} = cosa e cosX {y * } + sina e sinX { y Jt} cos«|> { y Jt} ; 


+ IF cos X ( sinX- cos(|). 


tana. 


\ijk} 


sml { iik } sm<\> u 


my 


r ik cos2 x j + r ; -,- cos2 h + r kj cos2 x i 


= 0. 


As reference 2 shows, this quadratic equation can be 
reduced to the following mutually exclusive solutions: 


Equation (15) has two solutions; however, unlike the 
solution for angle of attack, these solutions are not 
necessarily mutually exclusive. In general, selection of 
the correct root depends upon the port arrangement used 
to determine the angle of sideslip. For most practical 
arrangements, however, selection of the root of 
equation (15) whose magnitude is closest to zero will 


for |a e | < 45°, a e = ^tan , 

(13(a)) 

give the correct result. Reference 2 provides a detailed 
discussion of the sideslip root selection criteria. 

for|a e |>45°, a e = ^n-tan [g]j. 

(13(b)) 

Mach Number. Static Pressure, and Impact Pressure 
Estimator 


where the parameters A and B are defined as follows: 

.2 2 2 
A = F (/ r sin . X i + r . sin X k + F, siiF 

B = II* cosck sinX, cosX ■ + T ji cos^sin^cos^ (14) 

+ T kj cos(|) ; sinA cosL . 

Angle-of-Sideslip Estimator 

When the effective angle of attack has been estimated, 
then the effective angle of sideslip may be evaluated 
using any combinations of the available ports except the 
set where all three ports lie on the vertical meridian. The 
result is a quadratic equation in tan P , 


When the values of a g and \i ( , have been determined, 
the incidence angles at all of the ports can be evaluated, 
and only e, and q c remain unknown in the pressure 
model (eq. (10)). As described earlier, e typically is a 
function of Mach number, and hence is implicitly 
related to both p m and q ( .. As a result, solutions for p m 
and q c must be iteratively extracted at each time step. 

•J 

The resulting estimator" takes the form 


Pi 

Pi 


Vc 

(m + 1) 

Poo 

= 1- 


M (m) Q M (m) l M (m) Q 


A’ tan“P + 2B’ tanP + C' = 0, 


(15) 


LPnJ 


, (16) 


where 


where 


6 

American Institute of Aeronautics and Astronautics 



M 


(»?) 


( 2 A 

COS 0J 
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COS 0 ;j 
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+ e sin 0„ 
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. 0 

0 

q 2 
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. 0 

0 

0 

q 3 ■ 

. 0 

0 

0 

o . 
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In equation (16), the superscript (m) refers to the 
result of the w th iteration and n indicates the number of 
ports. For the first iteration, the value of e <m ' ) is set 
based on the final Mach number from the previous time 
step. The q ■ terms are weights that have a nominal value 
of 1.0. Setting the value of q ( to 0 weights the / th 
pressure reading out of the algorithm. Reference 2 
provides a detailed derivation of equation (16). 

When equation (16) is solved for q c and , Mach 
number can be computed using normal one-dimensional 
fluid mechanics relationships. For subsonic conditions, 
Mach number can be directly calculated using isentropic 
flow laws in which, assuming that the ratio of specific 
heats (y) is 1.4, 

q r 9 3.5 

- = [1.0 + 0.2M"] -1. (17) 

Po o 


For supersonic conditions, the solution is computed 
using normal-shock wave relationships. 14 For y = 1.4, 



166.92 M„ 


M„ 


7 M " - 1 


2.5 


- 1. 


(18) 


Equation (18) is easily solved using a Taylor's series 
expansion and a reversion of series to compute Mach 
number: 15 


jE r i (w ‘) 

-;W = 1.839371 (191 

where the coefficients of the summation are r = [ 1.42857, 
-0.357143, -0.0625, -0.025, -0.012617, -0.00715, 
-0.0043458, 0, 0, -0.0087725], 

For hypersonic Mach numbers, a very significant 
temperature rise (greater than 1000 °F) occurs across the 
bow shock wave and the ratio of specific heats can no 
longer be assumed to be constant. Therefore, some error 
is expected using equation (18) in the Mach number 
calculation. Fortunately, because pressure depends 
primarily upon the mechanical aspects of the flow, 
influences caused by high-temperature gas properties 
are secondary. Numerical analyses of equation (18) have 
shown that the error introduced is less than 0.2 percent 
for Mach numbers less than 4.0. 16 This error is 
considered acceptable for real-time applications for 
most aerospace vehicle designs. 

Flush Airdata Sensing Pressure Model 
Calibration Procedure 

The FPM described by equations (3) and (10) has 
three parameters that must be determined by empirical 
calibrations: 5a, 5(3, and e. These calibrations allow 
the FPM to be “shaped” to characterize pressure 
distributions and flow deflections for various forebody 
shapes. Given a set of reference conditions (obtained 
from wind-tunnel or flight data) that includes the surface 
pressure distribution, the true angles of attack and 
sideslip, the free-stream Mach number, and dynamic or 
static pressure, the model parameters may be estimated 
using the following algorithm: 

1. Equations (13) and (15) are used to estimate a g 
and |3 ( from the surface pressure distribution. The 
flow angle calibrations are then given by 
equations (4) and (5). 

2. Given a g and |3 ; . the surface incidence angles are 
evaluated for each of the pressure ports using 
equation (3). 
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3. The result of step 2 is substituted into 
equation (10) and rearranged: 

C p (0)-cos 2 e = sin 2 0 ■ e(M oo ,a e ,P ff ), (20) 

where 

ve) . *!>!£=. 

p c q c 

4. Equation (20) is evaluated and collecting for all of 
the n pressure ports yields 


Using the above procedures, the model calibrations 
for a new configuration can be systematically developed 
from the reference data. Typically, systematic trends in 
the calibration parameters are identified by plotting the 
estimated calibration parameters as a function of 
reference variables and visually inspecting the results. 
Care should be taken to ensure that adequate reference 
data are generated in order to cover the desired flight 
envelope. Also, calibrations with respect to different 
geometrical configurations for the vehicle (for example, 
gear up or down and control surface deflections) may be 
required. Detailed calibration examples for various 
forebody shapes are given in the “Results” section. 


C pc (0 1 )-cos 2 0 1 


sm 0, 


Blunt-Bodv Pressure Data Sources 


C pc (0„)- cos 0„ 


• -a 
sm 0, 


PJ- (21) 


5. Equation (21) can be directly solved 17 to give a 
least-squares estimate of e: 


X sin2 0;[c pc (0,-)-cos 2 0 ( J 

e (M 00 ,a g ,P g ) = ---. (22) 

X sii|l0 ( 

i = 1 

6. The estimate of £ derived from equation (22) 
gives the minimum fit error with respect to the 
pressure distribution. The estimate, however, may 
not give the minimum error with respect to the 
reference Mach number. If a refinement is desired 
to minimize the fit error with respect to Mach 
number, then a “steepest-descent” algorithm can 
be used. For M,. ef - M 0) > SM required , 

-i ( m ) 

£ m + 1 = £ (w ° + 3 4-[M,, f -M (n,) ]. (23) 

dM lel 


Three calibration parameters must be evaluated for 
any FADS system: the shape and compressibility 
parameter, e; the angle-of-attack flow correction angle, 
8a; and the angle-of-sideslip flow correction angle, Sp. 
As mentioned in the introduction, a general set of 
guidelines regarding the structure of the calibration 
curves for various FADS systems has never been 
developed and reported. This paper attempts to develop 
such a set of rules by deriving calibrations for a variety 
of shapes. The calibrations are derived from three 
distinct data sources: analytical solutions for simple 
geometries, wind-tunnel data, and flight data. The 
following subsections detail the geometrical 
configurations analyzed and the data sources. Results 
obtained by applying the FADS calibration procedures 
to the data are presented in the “Results” section. 

Analytical Solutions for Simple Geometries 

To better understand the variations in the FADS 
calibration caused by configuration changes, an 
evaluation of exact, analytical solutions was desired. 
Two fixed-geometry incompressible solutions already 
exist that exactly fit into the FPM. Incompressible, 
uniform crossflow over a cylinder has a pressure 
solution given by: 


In equation (23), (3e^ n0 )/3M is the sensitivity 
derivative of the calibration parameter with respect 
to Mach number, and M ,ml is the Mach number 
solution at the rn til iteration. The derivative 
(3 e^” 0 )/3M is numerically determined by 
perturbing Mach number and determining the 
change in e~" l \ 


_ P( Q )-Po o , . . 2 a 

C (0) = - = 1 - 4sm 0 


cos' 0-3 sin' 0. 


(24) 


Incompressible, uniform flow over a sphere has a 
pressure distribution described by: 
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„ P( Q )-P~ , 9 . 2 n 

C (0) = - = l--sm 0 

p <7 C 4 

= cos'O - ^sitCO. 

4 


(25) 


The FPM (eq. (10)), in which e is equal to -3 for a 
cylinder and -1.25 for a sphere, exactly models the 
incompressible flow solutions. 

A general, analytical, potential-flow solution 
algorithm 19 was developed for the three-dimensional 
ellipsoid shown in figure 2. This algorithm calculates 
the pressure coefficient at any point on the surface of the 
ellipsoid, and the ellipsoid can be arbitrarily inclined to 
the free stream. The ellipsoid fineness ratio, f, is defined 
as the ratio of the major to the minor ellipsoid axes; the 



major axis is also the reference axis. Analysis must be 
limited to conditions where flow separation is not 
expected because the incompressible solution will not 
predict it. 

Because the FPM was adapted from a spherical 
model, the ellipsoid provides a strenuous test of the 
FPM application to nonspherical shapes. This test is 
necessary if the FPM is to be generally applicable. The 
analytical ellipsoid pressure model allows the evaluation 
of the FADS calibration parameters with varying 
forebody shape, angle of attack, and angle of sideslip. 

Wind-Tunnel Data 

The potential flow solutions described in the previous 
section clearly do not include the effects of viscosity, 
Reynolds number, or compressibility. To understand 
these effects in greater detail, data from a series of wind- 
tunnel tests were analyzed. These data cover a wide 
range of Mach numbers, incidence angles, and forebody 
shapes. The wind-tunnel data also were used to 
demonstrate the sensitivity of the FADS calibration to 
the port locations. The test facilities, model shapes, test 
conditions, and wind-tunnel instrumentation are 
described hi detail in the following subsections. Table 1 
shows the range of test conditions analyzed for each of 
the wind-tunnel models. 

Forebodv Wind-Tunnel Tests 

A series of wind-tunnel tests, designed to test the 
adaptability of the FPM and the accompanying 
calibration procedures, was conducted in the National 
Technical Systems (Saugus, California) 4 ft-by-4 ft 


Table 1. Wind-tunnel model test conditions. 


Configuration 

Tunnel 

Mach number 

Low 

range 

High 

a range, 
deg 

Low High 

(3 range, 
deg 

Low High 

8-deg cone 

NTS 

0.30 

4.50 

-5 

30 

0 

0 

X-33 

forebody 

NTS 

0.30 

4.50 

-5 

30 

-5 

15 

Rankine 

half-body 

NTS 

0.30 

4.50 

-5 

30 

0 

0 

X-33 

16-Foot 

0.25 

1.20 

-10 

24 

-10 

16 

X-33 

UPWT 

1.60 

4.50 

-10 

24 

-10 

10 

X-34 

14- by 22-Foot 

0.30 

0.30 

-4 

20 

-8 

8 

X-38 

LMVS 

0.40 

0.95 

-5 

24 

-4 

4 
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blow-down-to-atmosphere wind tunnel. The facility 
operates at Mach numbers ranging from 0.25 to 5.00. 
High-pressure ah is fed through a variable-geometry 
nozzle to achieve the desired test conditions. The walls, 
floor, and ceiling in the transonic test section are 
22-percent porous to minimize reflected shock waves. 

This series of wind-tunnel tests was conducted with 
three models (fig. 3). The first model was a spherically 
blunted 8-deg nosecone. This shape is similar to several 
single-stage-to-orbit launch vehicle designs and tighter 
forebody shapes. The second model shape was achieved 
by rotating a two-dimensional Rankine half-body 
along its central axis. The result was an ellipse-like 
forebody that approaches a cylindrical shape at infinite 
distance. The third shape accurately modeled the first 
20 percent of body length of the X-33 single-stage-to- 
orbit vehicle (Lockheed Martin Skunk Works, Palmdale, 
California) at a scale of approximately 4.2 percent. A 
smooth boat-tail afterbody was added to the X-33 
forebody to move the base separation away from the 
nosecap region. This model represents 
nonaxisymmetric, lifting-body type configurations. 

Each of the models consisted of three parts: the 
nosecap, the body, and a common mounting adapter. 
The nosecap was made of stainless steel to minimize 
any damage during the wind-tunnel blows. Figure 4 




**The Rankine half-body conies from the inviscid, incompressible 
solution for a point source in uniform flow. 


shows the port layout and numbering system for the 
nosecap. Each of the FADS pressure ports were located 
and drilled normal to the surface. Stainless steel tubing 
(0.031 in. inside diameter) connected the external model 
surface to the lubber tubing that was plumbed to the 
pressure transducer. The model bodies were made of 
aluminum. A stainless-steel mounting adapter served as 
the structure that held the nosecap and body, mounted 
the completed model to the wind-tunnel sting, and 
housed the model instrumentation. 

Two types of instrumentation were mounted in the 
model adapter. First, a 32-port electronic scanning 
pressure (ESP) module (Pressure Systems Incorporated, 
Hampton, Virginia) was located close to the nosecap to 
measure differential pressure (±15 lbf/in') between the 
nose ports and atmospheric pressure. The second 
internal instrumentation was a digital inclinometer, used 
for measuring the model incidence angle during the 
wind-tunnel runs. Because no balance was used in the 
test, the inclinometer was used to account for sting 
bending under loaded conditions. When the model was 
upright, the inclinometer measured angle of attack; 
when the model was rolled 90°, the inclinometer was 
reinstalled to measure angle of sideslip. 

X-33 Wind-Tunnel Tests 

FADS wind-tunnel calibration data were obtained for 
a 2-percent model of the X-33 vehicle (designation 
604B002C) in the 16-Foot Transonic Tunnel 20 and the 
Unitary Plan Wind Tunnel at NASA Langley. The 
16-Foot Transonic Tunnel is a closed-circuit, single- 
return, continuous-flow, atmospheric tunnel that 


990229 

Figure 4. Port arrangement of wind tunnel-test models. 


Port 

A., deg 

<t>, deg 

1 

0 

0 

2 

20 

0 

3 

20 

90 

4 

20 

180 

5 

20 

270 

6 

30 

0 

7 

30 

90 

8 

30 

180 

9 

30 

270 

10 

45 

0 

11 

45 

90 

12 

45 

180 

13 

45 

270 

14 

60 

0 

15 

60 

90 

16 

60 

180 

17 

60 

270 
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operates at Mach numbers ranging from 0.2 to 1.3. The 
tunnel has an octagonal. 16-ft slotted test section. Some 
run limitations are sometimes imposed near transonic 
speeds because of the lack of porous walls to absorb 
reflected shock waves. The small size of the X-33 
model, however, allowed testing at transonic Mach 
numbers without reflected shock-wave interference. The 
Unitary Plan Wind Tumiel is a closed-circuit, 
continuous-flow, variable-density, supersonic wind 
tunnel with two 4 ft—b\'—4 ft test sections. One test 
section operates at a Mach number range from 1.5 to 
2.9, and the other operates at a range from Mach 2.3 to 
approximately Mach 4.5. The Reynolds number was 
nearly constant at 2xl0 6 /ft. Results of the X-33 

•y 

FADS wind-tunnel tests were also previously reported." 

The 2-percent X-33 model was built to measure total 
vehicle forces and moments and nosecap surface 
pressures. The wind-tunnel model was machined from a 
solid piece of aluminum, and the nosecone was made of 
stainless steel to minimize damage from tunnel 
contaminants. Tests were conducted with and without a 
grit ring installed aft of the FADS ports, but little 
variation existed in the pressure distribution. 

The nosecap surface pressures (fig. 5) were sensed 
using an onboard ESP module that produced a time- 
multiplexed analog output with 32 channels. The ESP 
sensor was a ±10-lbf/in 2 differential module with a 
manufacturer’s accuracy of better than ±0.1 percent of 
full-scale reading. The time-multiplexed analog outputs 
from the ESP module were tagged and sampled using a 
16-bit analog-to-digital conversion system. The 
2-percent X-33 model FADS pressure ports were drilled 
normal to the surface with a diameter of 0.02 in. 

X-34 Wind-Tunnel Tests 

The FADS wind-tunnel calibration data were obtained 
for a 10-percent scale model of the X-34 vehicle 
(Orbital Sciences Corporation, Dulles, Virginia) in the 
14- by 22-Foot Subsonic Tunnel at NASA Langley. The 
14- by 22-Foot Subsonic Tumiel"" features a closed 
circuit and operates at atmospheric pressure. The X-34 
tests were performed at a nominal Mach number of 0.3 

-y 

and a dynamic pressure of 125 lbf/ft". The free-stream 
angle of attack ranged from -4° to 20°; the angle of 
sideslip ranged from -8° to 8°. 

The model was made of a carbon fiber-fiberglass 
composite and was sting-mounted at the base. Figure 5 
shows the locations of the ports. The X-34 model FADS 
orifices were connected to a 32-port ESP module. The 

•y 

ESP sensor was a ±2.5-lbf/in differential module with a 
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Figure 5. X-33, X-34, and X-38 port arrangements. 


manufacturer’s accuracy of better than ±0.1 percent of 
full-scale reading. The time-multiplexed analog outputs 
from the ESP module were tagged and sampled using a 
16-bit analog-to-digital conversion system. The X-34 
model FADS pressure ports with a diameter of 0.06 in. 
were drilled normal to the surface. ESP zero-shift 
calibrations were performed when the tumiel 
temperature changed by more than 4.5 °F and after each 
model changeover. 

X-38 Wind-Tunnel Tests 


The X-38 vehicle (designated VI32) is similar to the 
X-24A lifting-body design. The X-38 deploys and lands 
under a parafoil rather than on a runway. Modifications 
to the X-24A design include removing the center fin to 
allow for parafoil storage, redesigning the aft structure 
to improve the base pressure, and replacing the booster 
rockets with a drogue chute. 


11 

American Institute of Aeronautics and Astronautics 






The nine-hole FADS system installed on the nosecap 
(fig. 5) was calibrated at the Lockheed Martin Vought 
Systems (LMVS) transonic wind tunnel in Dallas, 
Texas. The 4 ft-by-4 ft blow-down-to-atmosphere wind 
tunnel operates at Mach numbers ranging from 0.25 to 
5.00. High-pressure air is fed through a variable- 
geometry nozzle to achieve the desired test conditions. 
A 5.2-percent scale model was used to represent the 
23-ft full-scale X-38 vehicle. The Mach range varied 
from 0.40 to 0.95. The angle of attack varied from -5° 
to 24°, and the angle of sideslip varied from -4° to 4°. 

Miscellaneous Wind-Tunnel Test Data 

A search for literature containing experimental 
pressure distribution data on generic shapes yielded 
information that contributes to the goals of the current 
study. Data were found for a sphere; spherically 
blunted 9-deg, 12.5-deg, 18-deg, and 45-deg 

nosecones; 24, 25, 26 and a spherically capped cylinder. 
The 12.5-deg, 18-deg, and 45-deg nosecone data were 
only available at hypersonic Mach numbers (greater 
than Mach 6), because these shapes were being 
considered for reentry capsule shapes. Most reentry 
capsule experiments were exclusively conducted at 0° 
angle of attack. 

Flight Test Data 

Data obtained from two full-scale FADS flight tests 
are presented in the following subsections. These flight 
data were obtained from a real-time FADS 
demonstration flown on the NASA Dry den F/A-18 
SR A 10 and from transonic FADS tests on an F-14 

o 

aircraft. The flight data results are compared with 
analytical and wind-tunnel analyses in the “Results” 
section. 

F/A-18 Systems Research Aircraft 

Flush Airdata Sensing Flight Tests 

The SRA flight tests demonstrated a prototype real¬ 
time FADS system. The test envelope included Mach 
numbers as high as 1.6 and angles of attack greater than 
45°. The airdata were sensed by 11 orifices located on a 
composite nosecap molded into the structure of the 
aircraft radome and faired smooth to the surface. 
Figure 6 shows the locations of the pressure ports on the 
nosecap. 

The SRA FADS pressures were sensed by 11 
miniaturized, digital, absolute-pressure transducers. 


Port 

K deg 

<t>> deg 

1 

0 

0 

2 

40 

0 

3 

40 

180 

4 

55 

0 

5 

55 

90 

6 

55 

180 

7 

55 

270 

8 

60 

45 

9 

60 

135 

10 

60 

225 

11 

60 

315 


Port 

/l, deg 

<(>, deg 

1 

60 

180 

2 

40 

180 

3 

20 

180 

4 

0 

0 

5 

20 

0 

6 

40 

0 

7 

60 

0 

8 

60 

90 

9 

30 

90 

10 

30 

270 

11 

60 

270 


F-18/SRA 


\l 10 

-fr-T 

t 2 . 

T 4 i*i 


F-14 

tl 

*2 

8 9 T 3 10 11 

*5 
• 6 

*7 
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Figure 6. F/A-18 SRA and F-14 port arrangements. 


Each absolute-pressure transducer incorporates a four- 
active-arm strain-gage bridge for data sensing; internal 
signal conditioning with a 20-Hz antialiasing low-pass 
filter and output signal amplification; and a 20-bit 
analog-to-digital conversion. The transducers have a 
repeatability that exceeds 0.01 percent of full scale and 
a measurement range from 1.50 to 40.00 lbf/in 
atmosphere. Data collection and algorithm 
computations are performed onboard in real time by two 
commercially available single-board computers inserted 
in a ruggedized flight chassis. The raw pressure 
measurements and onboard airdata computations were 
telemetered to a ground station where they were 
recorded for postflight analyses. 

The SRA FADS system was calibrated by comparing 
the onboard calculations with flight-derived reference 
airdata values and adjusting the calibration coefficients 
using the previously described calibration procedures. 
The flight reference data were generated by merging 
complementary information from multiple data sources, 
including the onboard inertial navigation system 
attitudes, rates, and accelerations; radar-tracking 
velocity and position data; and raw insonde weather 
balloon data. 
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F-14 Flush Airdata Sensing Flight Tests 

The F-14 FADS flight tests were performed at NASA 
Dryden on an aircraft with special equipment 
modifications that were part of a program demonstration 
for new flight control techniques.” The equipment 
modifications included a pair of deployable, canard-like 
surfaces on the fuselage forebody for spin recovery and 
a gun fairing in the left side of the fuselage forebody. 
Although the canards were closed for the FADS tests, 
they may have some effect on the upstream flow. The 
geometric asymmetry caused by the gun fairing may 
have also influenced some of the FADS pressure 
measurements. 

For the F-14 tests, 11 flush orifices (0.031 in. 
diameter) were installed on a spherically blunted 
nosecap in a cruciform pattern. The nosecap blended 
with the nose of the aircraft. Orifices along the vertical 
axis were placed at angular increments of 20°; those 
along the horizontal axis were placed at increments of 
30°. Additional flush-mounted orifices were installed on 
the nose section aft of the nosecap. These 0.125 in- 
diameter orifices were located in 2 rows around the 
fuselage: a row of 4 orifices evenly spaced 24.0 in. aft of 
the nose apex and another row of 12 orifices 50.0 in. aft 
of the nose apex. The nose-section orifices were 
installed to furnish a suitable static-pressure source to be 
used in conjunction with the orifices on the nosecap. 
Figure 6 shows the locations of the FADS ports. Details 
of the F-14 instrumentation system are reported in 
reference 8. 

For the F-14 FADS tests, data were obtained for Mach 
numbers between 0.60 and 1.60, for angles of attack to a 
maximum of 26.0°, and for angles of sideslip to a 
maximum of 11.0°. The flight Reynolds numbers varied 
between 1.0xl0 6 /ft and 5.5xl0 6 /ft. High flow 
angles at high transonic speeds were limited. 

Results 

The suitability of various vehicle shapes for a 
FADS system is discussed in this section. First, 
analytical solutions for the incompressible 
pressure field over a sphere and varying fineness-ratio 
ellipsoids are used to validate the suitability of the FPM 
and calculate the appropriate FADS calibration 
parameters (e, 8a, and SP) at incompressible flow 
conditions. Varying the ellipsoid fineness ratio reveals 
the behavior of the FPM calibration for varying 
forebody shape; varying the free-stream angle of attack 
reveals the behavior of the calibration to changes in the 
vehicle orientation to the free stream. 


Second, using wind-tunnel and flight data, the 
variation in the calibrations is extended to a wide variety 
of blunt-forebody shapes. The effect of shape and 
compressibility is assessed at Mach numbers to a 
maximum of 5.0. 

Last, the sensitivity of the FADS calibration to the 
port locations is addressed. The modeling suitability of 
the FPM is illustrated by plotting the impact pressure 
coefficient, C , as a function of the incidence angle. 

rC 

The C is found by rearranging the FPM (eq. (10)): 

r c 

W(0) - p 9 9 

C (0) = - — = cos'0 + esin'0. (26) 

p c q c 

Ellipsoid Analytical Data 

The ellipsoid pressure calculations provide a 
comprehensive data package to test the validity of the 
FPM. Although the data represent only incompressible 
flow conditions, variations in the FADS calibration 
caused by forebody shape and angle of attack were 
investigated. 

Validation of the Flush Airdata Seasing 

Pressure Model 

The ability to model the pressure distribution over the 
forward portion of a blunt body is essential to ensure 
general applicability of the FPM to a variety of vehicles. 
The analytical ellipsoid solution 19 was used to generate 
forebody pressures for ellipsoids ranging from a sphere 
(f=l) to a high-fineness-ratio ellipsoid (f=20). Figure 7 
shows the calculated impact-pressure coefficients and 
the calibrated FPM as a function of the flow incidence 
angle, 0, for varying ellipsoid fineness ratios at 0° angle 
of attack. In all cases, a calibration parameter, e, exists 
that, when used in equation (10), identically matches the 
analytical pressure distribution. As expected (eq. (25)), 
the value of e for a sphere is equal to -1.25. As the 
fineness ratio becomes large, the value of £ 
approaches 0. 

As the angle of attack increases to 30°, no degradation 
in the pressure modeling capability of the FPM is found. 
Figure 8 shows that the least-squares estimate of £ 
results in a perfect match of the data. Note that the 
calibration for the sphere is unchanged with increasing 
angle of attack; however, the nonspherical shapes 
exhibit a reduced sensitivity of the calibration to 
fineness ratio at angles of attack greater than 10°. 
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Calibration of the Shape and Compressibility 
Parameter 


-o— 8 =-.014, f = 20.00 
8 = -.122,f = 5.00 
-O- 8 = - .464, f = 2.00 
-*- 8 = - .700, f = 1.50 
8 = - .908, f = 1.25 



Figure 7. Comparison of the FPM and the analytical 
impact-pressure coefficient as a function of flow 
incidence angle, ellipsoid, a = 0°. 


-o- 8 = - .747, f = 20.00 
8 = - .824, f = 2.00 
-o- 8 =-.933,1 = 1.50 
8 = -1.050,1 = 1.25 



Figure 8. Comparison of the FPM and the analytical 
impact-pressure coefficient as a function of flow 
incidence angle, ellipsoid, a = 30°. 


Because the FADS solution algorithm computes 
effective angles of attack and sideslip first, the e 
variation caused by angle of attack should be 
documented as a function of effective angle of attack 
rather than true angle of attack. This method also 
reduces the nonlinearity in the calibration model. 

Figure 9 shows a plot of e as a function of effective 
angle of attack for varying fineness-ratio ellipsoids. The 
plot confirms that e is equal to -1.25 for a sphere (f= 1). 
Because the flow field of the sphere is independent of 
orientation, £ does not vary with changes in angle of 
attack. For small increases in the fineness ratio (for 
example, f=1.25), £ significantly increases and the 
dependence of £ on the effective angle of attack 
becomes noticeable. Increasing fineness ratio results in 
greater concavity of the E-as-a-function-of-og, curve. 
As the fineness ratio of the ellipsoid increases, causing 
its shape to approach a slender rod, the value of £ 
approaches 0 at 0° angle of attack. 


-1 = 1.00 1 = 2.50 

-1 = 1.25-1 = 3.00 

-1 = 1.50-1 = 4.00 

-1 = 1.75 1 = 8.00 



Figure 9. FPM calibration variation with effective angle 
of attack for varying fineness-ratio ellipsoids. 


Evaluating the FPM (eq. (10)) for £ = 0 at X = 0° 
and K = 90° yields total and static pressure, 
respectively. These results match the expected 
measurements of total pressure at the tip (X = 0°) and 
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static pressure along the side (X = 90°) of a thin rod in 
subsonic flow. The shape of the e = e(a g ) curve can 
be approximated by a polynomial equation, allowing 
simple calibration storage in a real-time airdata system. 

Calibration of the Upwash Correction 

Equation (4) requires an upwash correction, 8a, to 
compute true angle of attack from effective angle of 
attack. Figure 10 shows a plot of the upwash correction 
for varying fineness-ratio ellipsoids in incompressible 
flow. As expected for a sphere, the upwash correction is 
0. The upwash correction is also 0 for all axisymmetric 
ellipsoids at 0° angle of attack. 

As the angle of attack increases, the upwash 
correction becomes larger, peaking between 45° and 53° 
effective angle of attack (depending on fineness ratio). 
The upwash correction decreases beyond the peak angle 
of attack. As fineness ratio increases, the upwash 
increases as well, reaching an apparent upper limit for 
an ellipsoid with a fineness ratio greater than 15. As was 
the case with e, the upwash parameter can be 
approximated with a polynomial in effective angle of 
attack for a given Mach number. 

The sensitivity of the upwash correction to variations 
in angle of sideslip was investigated for ellipsoids with 


-f = 1.00 f = 2.50 

-f = 1.25-f = 3.00 

-f = 1.50-f = 4.00 

-f = 1.75 1 = 8.00 


- 1 = 2.00 . 1 = 20.00 



a e , deg 

990235 

Figure 10. Variation of upwash correction with effective 
angle of attack for varying fineness-ratio ellipsoids. 


varying fineness ratios. Figure 11 shows the angle-of- 
attack error introduced, using the calibration generated 
for 0° angle of sideslip, at a nonzero angle of sideslip at 
30° angle of attack. In all cases, the angle-of-attack error 
is negligible. The error can be larger for 
nonaxisymmetric bodies, as is shown in the next 
section. 

- 1 = 1.00 

-1 = 1.25 

-1 = 1.50 

- 1 = 2.00 

-1 = 3.00 


1 = 5.00 
1 = 10.00 
1 = 20.00 



Figure 11. Angle-of-attack error introduced by 
neglecting the sideslip effect on the e calibration. 


Calibration of the Sidewash Correction 

Because the ellipsoid is axisymmetric, the sidewash 
correction at 0° angle of attack is identical to the upwash 
correction at 0° angle of sideslip. The upwash correction 
was previously shown to be independent of sideslip for 
moderate sideslip angles. The sidewash correction, 
however, can be a significant function of angle of attack. 

Figure 12 shows the error obtained using the ellipsoid 
0° angle-of-attack sidewash calibration at 30° angle of 
attack. As the ellipsoid deviates from a sphere, the 
sideslip error increases, reaching 25 percent for a 
fineness ratio of 3.0 at 5° angle of sideslip. As a result, 
most wind-tunnel calibration efforts will require both 
angle-of-attack and -sideslip sweeps to gather sufficient 
data. 
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Figure 12. Sideslip error introduced by neglecting 
angle-of-attack effect on the e calibration. 

Blunt-Bodv Wind-Tunnel and Flight Data 

A study of the FPM using the analytically calculated 
ellipsoid pressure data has revealed the general behavior 
of the calibration with varying vehicle shape and 
orientation to the free-stream flow field. The analytical 
model, however, lacks important flow features that most 
atmospheric vehicles encounter: Mach number 
(compressibility), Reynolds number, and flow 
separation effects. Data from a series of wind-tunnel 
experiments and from flight tests are used in the 
following section to evaluate the FADS model 
characteristics for a variety of vehicles at realistic flight 
conditions. 

Effect of Flow Separation 

The FPM is based on potential flow over a sphere. 
This model assumes that the fluid flow is always 
attached over the region of interest, primarily where the 
ports are located. Wind tunnel-derived pressure 
distributions on a sphere at two Reynolds numbers are 
used to study the breakdown of the FPM under 
separated flow conditions. 

Figure 13 shows pressure distributions measured on a 

no 

sphere at a Reynolds number based on a diameter of 
approximately 4.5 X 10 ; the potential flow solution for 
a sphere; and the calibrated FPM. At this Reynolds 
number, the flow stays attached to an incidence angle of 


O Rel. 23 
— FADS pressure 



Figure 13. Comparison of FPM and wind-tunnel 
impact-pressure coefficient on a sphere, 
Re^ = 4.2 to 4.6 x 10 5 . 

more than 110° and the calibrated FPM does an 
excellent job of modeling the pressure distribution. 

At lower Reynolds numbers (Re^ = 1.6 X 10 5 ), 
however, a laminar separation occurs at approximately 
60° incidence angle (fig. 14). Attempting to calibrate the 
FPM for the 0°-l 10° range of incidence angle results in 
a poor match. Restricting the incidence-angle range to 


O Rel. 23 

— FPM calibrated using ports 
in attached flow 

— FPM calibrated using ports 


in separated flow 



Figure 14. Comparison of the FPM and the wind- 
tunnel impact-pressure coefficient on a sphere, 
Re^ = 1.6 to 1.7 X 10 5 . 
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the ports in the attached flow region results in a good 
model of the pressure distribution. Using wind-tunnel 
data to identify the flight conditions at which individual 
ports have separated flow is generally conservative, 
because the full-scale flight data usually are at higher 
Reynolds numbers (delayed separation). 

Two methods exist of enforcing the restriction on 
using ports in separated flow regimes. One method is to 
locate ports so that separation does not occur in the 
flight envelope of interest. For the X-33 vehicle, the 
largest port cone angle on the leeward vertical meridian 
is 20°, which does not have separated flow to a 
maximum of 45° angle of attack. An alternate port 
arrangement method is to adjust the ports that are used 
in the FADS algorithm based on flight condition. This 
method is complex, possibly requiring multiple 
calibrations and data smoothing during transitions 
between port sets. 

Validation of the Flush Airdata Sensing 

Pressure Model 

Pressure distributions measured during the wind- 
tunnel tests of the 8-deg cone, X-33 forebody, and 
Rankine half-body were used to validate the pressure 
modeling capability of the FPM from Mach 0.40 to 
Mach 4.75 from -4° to 30° angle of attack. Because 
some of the ports on the nosecap were hr or near 
separation zones at high angles of attack, a realistic 
subset of the available ports was used in the calibration 
algorithm to estimate e. Experience on F-18 and X-33 
FADS systems 2, 3 ’ 10 and analysis of the wind-tunnel 
data were used to select a reasonable set of ports: 1,2,3, 
4, 5, and 10 (fig. 4). Other combinations of ports could 
be used; however, the simplified FADS solution 
algorithm described in this paper requests a minimum of 
five ports, three of which must be in a vertical line for 
angle-of-attack estimation. 

Figure 15(a) shows pressure distributions over the 
axisymmetric, spherically blunted 8-deg nosecone 
configuration and the calibrated FPM for Mach numbers 
of 0.4, 0.9, 1.2, and 3.5 at 0° angle of attack. In each 
case, the pressure distribution is well-modeled to greater 
than 45° incidence angle. A small difference in the 60° 
incidence-angle ports—which were not used to calibrate 
the shape—is noted. A challenging test for the model is 
to accurately model the high-angle-of-attack pressure 
distributions. Figure 15(b) shows the pressure 
distributions and FPM fit at 30° angle of attack for 
subsonic to high supersonic Mach numbers. The FPM 
accurately models the pressure distribution to a 
minimum of 60° flow incidence angle for all Mach 
numbers. 


Unlike the ellipsoid and the 8-deg cone, the X-33 
forebody model is not axisymmetric. Figures 16(a) and 
16(b) show the FPM fit of the measured X-33 forebody 
pressure distributions for Mach numbers of 0.4, 0.9, 1.2, 
and 3.5 at angles of attack of 0° and 30°, respectively. 
As with the 8-deg nosecone, the pressure distribution is 
well-modeled to 45° incidence angle, with some 
degradation occurring at 60°. At 30° angle of attack, the 
pressure distribution is well-modeled to greater than 60° 
incidence angle. No degradation in the FPM predictive 
capability was found for any of the Mach numbers or 
angles of attack tested. 

Although axisymmetric, the nose of the Rankine half¬ 
body does not have a spherical shape. The ellipsoid data 
previously shown confirmed that the FPM was 
adaptable to nonspherical shapes in incompressible flow 
fields. Figures 17(a) and 17(b) show the pressure 
distributions measured on the nosecap of the Rankine 
half-body wind-tunnel model for Mach numbers of 0.4, 
0.9, 1.2, and 3.5 at angles of attack of 0° and 30°, 
respectively. Because this model violates the spherical 
shape upon which the FADS model is based, some 
degradation in the predictive capability might be 
expected. This expectation is not the case, however, as 
the fit is shown to be as good as with either of the two 
spherically blunted nosecone models. These data 
confirm that the FADS pressure model can be 
successfully applied to nonspherical shapes from 
subsonic to supersonic Mach numbers, thereby 
extending the incompressible ellipsoid study findings to 
higher Mach numbers. 

Calibration of the Shape and Compressibility 

Parameter 

Figures 18(a)-(d) show the e calibrations for the 
spherically blunted 8-deg nosecone, Rankine half-body, 
X-33 vehicle, and X-33 forebody. The character of the 
calibration is similar for each forebody shape. At 
subsonic Mach numbers, e decreases with increasing 
effective angle of attack. For the axisymmetric 8-deg 
cone (fig. 18(a)) and Rankine half-body shapes 
(fig. 18(b)), the peak of the curve is at 0° angle of attack 
for subsonic and transonic Mach numbers. This peak 
location is expected because the flow behavior is 
symmetrical for positive and negative incidence angles 
on an axisymmetric shape. The X-33 vehicle, however, 
is not symmetric along the upper and lower surfaces aft 
of the spherically blunted nosecap. This asymmetry 
results in a peak e value at a nonzero angle of attack 
that is especially evident for the complete X-33 
configuration (fig. 18(c)). 
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(a) a = 0°. (b) a = 30°. 

Figure 15. Comparison of the FPM and the measured impact-pressure coefficient as a function of flow incidence 
angle, 8-deg cone. 




(b) a = 30°. 


Figure 16. Comparison of the FPM and the measured impact-pressure coefficient as a function of flow incidence 
angle, X-33 forebody. 
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(a) a = 0°. 


(b) a = 30°. 


Figure 17. Comparison of the FPM and the measured impact-pressure coefficient as a function of flow incidence 
angle, Rankine half-body. 


o M = 0.40 ► M = 1.50 
°Mr0.60 * M = 2.00 

o M = 0.90 4 M = 2.50 

v M = 0.95 » M = 3.50 

* M = 1.05 + M = 4.50 

M = 1.20 



(a) 8-deg cone calibration. 


o M = 0.40 v M = 0.95 * M = 2.00 

□ M = 0.60 * M = 1.05 » M = 2.50 

• M = 0.80 <M = 1.20 + M = 3.50 



(b) Rankine half-body calibration. 


Figure 18. FPM calibrations. 
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For Mach numbers less than 1.0, the angle of attack of 
the e peak increases with increasing Mach number. A 
comparison of figures 18(c) and 18(d) reveals 
differences between the X-33 and X-33-forebody 
calibrations at subsonic and transonic Mach numbers. 
Figure 18(e) shows a comparison of the complete and 
forebody-only X-33 configurations at Mach 0.60, 0.95, 
and 3.50. The calibration differences indicate that 
calibration of a blunt-body FADS system requires a 
complete configuration be tested to ensure that adequate 
data are obtained. 

Although the different forebody shapes have different 
values of e at low Mach numbers, all calibrations 
converge to a value near 0 at supersonic Mach numbers. 
As stated earlier, the FPM reduces to Newtonian impact 
theory as e approaches 0. The calibrations show that 
Newtonian impact theory is an accurate predictor of the 
surface pressures for these blunt bodies beginning at 
approximately Mach 2.0. The second-order nature of the 
£ calibration with angle of attack reduces such that at 
high supersonic Mach numbers, £ is nearly independent 
of angle of attack. 

Figure 19 shows a comparison of the e calibration for 
three (X-33, X-34, and X-38) FADS systems in 
development for low subsonic speeds. The X-33 data are 
front the wind-tunnel test of the complete X-33 
vehicle. 2 Both the X-33 and X-38 lifting bodies have 
camber that results in a peak in the £ curve at an 
effective angle of attack slightly greater than 0°. The 
X-34 plot, however, has an £ curve that peaks at an 
effective angle of attack of approximately 19°. Although 
the FADS ports are located symmetrically around the 
X-34 axisymmetric forebody, the forebody itself is 


rotated down 16.1° from the vehicle reference line, 
causing the shift in the £ peak. The most important 
feature of the £ calibration is that it is adaptable to a 
wide variety of realistic forebody shapes. 

To further study the Mach number effect on the 
calibrations, e m was defined as the peak of the curves 
of £ as a function of a g . Figure 20 shows e m as a 
function of Mach number for the 8-deg cone; Rankine 
half-body; X-33 forebody; X-33, X-34, X-38, F/A-18 
SR A, and F-14 vehicles; hemispherical-head cylinder; 
and 12.5-deg, 18-deg, and 45-deg spherically blunted 
cones. In all cases, the value of £ is a smooth, 
continuous, monotonic function of Mach number. This 
characteristic is an important difference from traditional 
pitot-static airdata system calibrations, which typically 
are characterized as having sharp changes in the 
position error calibration 1 in the transonic Mach region. 

An obvious outcast in figure 20 is the e m calibration 
for the F/A-18 SR A. The value of e is positive for low 
Mach numbers, rises slightly in the transonic region, 
and then decreases. If calibration data were achievable 
at high Mach numbers on the F/A-18 SR A, £ M would 
be expected to converge to zero (modified Newtonian 
theory). The likely reason for the difference in shape of 
the e calibration is that the region over which the FADS 
pressures operate on the SRA is in the first 2 in. of the 
forebody. The large afterbody of the F/A-18 SRA in 
relation to the forebody region causes higher pressure in 
the nosetip area, resulting in a higher e M . The result of 
the foregoing discussion is that traditional fighter 
forebody configurations may have different calibration 
shapes than blunt-body configurations (for example, 



a deg 
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-o- 8-deg cone 
-d- Rankine half-body 
-A- X-33 forebody 
—o— X-33 

-v- Hemispherical-head cylinder 
-ffl- 9-deg cone 
• X-34 


■ X-38 

▲ 12.5-deg cone 
♦ 18-deg cone 
▼ 45-deg cone 
N Sphere 

-F/A-18 SRA 

-F-14 



Figure 19. Comparison of X-33, X-34, and X-38 Figure 20. Variation of FPM calibration with Mach 

FPM calibrations, M ~ 0.3. number. 
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X-33, X-38, and Space Shuttle configurations); 
however, the FPM still yields results with equivalent 
levels of accuracy. 

Figure 21 shows the variation of e with angle of 
sideslip for the X-33 vehicle. At subsonic Mach 
numbers, e slightly depends on (3. Over the angle-of- 
sideslip range of interest (less than 5° for the X-33 
vehicle), the adjustment caused by angle of sideslip is 
small. The angle-of-sideslip effect is negligible at 
transonic and (as expected) supersonic Mach numbers. 

° Moo = 0.25 >M oo = 1.05 

oM oo = 0.60 *11/1^ = 1.10 

v Moo = 0.80 * Moo = 1-60 

a Moo = 0.90 « Moo = 2.00 
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Figure 21. X-33 FPM calibration variation with sideslip. 


Calibration of Upwash Correction 

An upwash calibration parameter, 5 a, is needed to 
calculate true angle of attack from the effective (locally 
measured) angle of attack, a e , according to 
equation (4). Figure 22 shows the upwash calibration 
for the three forebody models. All three blunt-body 
configurations have upwash calibrations that increase 
with increasing effective angle of attack to a peak value 
between 25° and 40° effective angle of attack at 
subsonic and transonic Mach numbers. Although data 
are limited, the upwash correction appears to decrease at 
high effective angles of attack, as was seen with the 
ellipsoid shapes (fig. 10). 


-M = 0.40 

-M = 0.60 

-M = 0.90 

-M = 0.95 

-M = 1.20 

.M = 2.50 



(a) 8-deg cone calibration. 
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-M = 0.80 
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(b) Rankine half-body calibration. 


Figure 22. Upwash calibrations. 


22 

American Institute of Aeronautics and Astronautics 



























































-M = 0.40 

-M = 0.60 

-M = 0.80 

-M = o.90 

-M = 0.95 

.M = 1.20 

-M = 1.80 

- M _ 2,50 



a_, deg 

990256 


(c) X-33 forebody calibration. 

Figure 22. Concluded. 

Figure 23 shows the upwash calibration for the 
complete X-33 configuration. The calibration is a 
smooth, polynomial function of effective angle of attack 
that decreases with increasing Mach number. The plots 

o Moo = 0.25 > Moo = 1.05 
o Moo = 0.60 * Moo = 1.20 
vM oo = 0.80 * Moo = 1.60 
a Moo = 0.90 -Moo = 2.00 



Figure 23. X-33 upwash calibration. 


show that the upwash correction reduces to 
approximately 0.25° at Mach 1.6 and greater. This 
reduction is expected at supersonic Mach numbers 
because the flow field is not disturbed by the model until 
the shock wave is encountered. The actual correction 
that will be used by the X-33 vehicle will normally be 
less than 3° based on the current Mach number and 
angle-of-attack trajectories. 

An important operational consideration is whether the 
FADS calibration can be easily implemented into a real¬ 
time estimator. Because the upwash calibration is a 
smooth, monotonic function of effective angle of attack, 
it can be represented as a simple polynomial or table 
lookup. 

A comparison between the X-33 and X-33 forebody 
upwash corrections (figs. 22(c) and 23) reveal 
differences of 1° in some cases. Again, this result 
indicates that the complete vehicle configuration must 
be used during blunt-body FADS calibration wind- 
tunnel tests. 

Figure 24 shows the upwash correction for the X-33, 
X-34, and X-38 forebodies at low subsonic Mach 
numbers. The slope of the upwash correction indicates 
that the X-38 model generated the largest upwash of the 
set. At 20° true angle of attack (approximately 37° 
effective angle of attack), the X-38 upwash correction is 
a large 17°. In contrast, the X-33 upwash correction at 
20° true angle of attack (approximately 27° effective 
angle of attack), is only approximately 7°. The larger 
X-38 upwash is caused by the larger thickness ratio of 
the fuselage. 
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Figure 24. Comparison of X-33, X-34, and X-38 
upwash calibrations, M = 0.3. 
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Calibration of Sidewash Correction 

For the axisymmetric shapes, 8(3 at 0° angle of attack 
is exactly the same as 8a at 0° angle of sideslip. For 
nonaxisymmetric shapes such as the X-33 vehicle, the 
sidewash calibration will obviously differ. Figure 25 
shows the sidewash calibration, 8(3, for the complete 
X-33 configuration. The sidewash calibration is similar 
in character to the X-33 upwash calibration (fig. 23), 
although the magnitude is approximately one-half of the 
upwash calibration at subsonic Mach numbers. The 
upwash correction should be larger than the sidewash 
correction, because the vehicle develops more lift with 
angle-of-attack changes than it does sideforce with 
angle-of-sideslip changes. As with the upwash 
calibration, the important feature is that the sidewash 
calibration is a smooth, monotonic function that is 
easily represented by a polynomial equation. 


0 1^ = 0.25 > Moo = 1.05 
o Moo = 0.60 * Moo = 1.10 
v Moo = 0.80 * Moo = 1.60 
a Moo = 0.90 « Moo = 2.00 



Figure 25. X-33 sidewash calibration. 


For the X-33 FADS calibration, 2 the error in 
neglecting the angle-of-attack effect on the angle of 
sideslip was less than 0.5° (la) from Mach 0.2 to 
Mach 4.5. Tims, the dependence of the sidewash 
calibration on angle of attack is configuration- 
dependent. A configuration that has a significant angle- 
of-attack effect on the sidewash will require 
considerably more wind-tunnel runs to calibrate than 
one that does not. 

Port Sensitivity 


but independent of the port locations. This idea came 
from early tests of the F/A-18 HARV, which showed 
only minor changes in the calibration with changes in 
the ports used in the model. Current wind-tunnel tests of 
the X-33 vehicle, X-33 forebody, 8-deg cone, and 
Rankine half-body, however, have shown the upwash 
calibration is somewhat dependent on the port layout for 
blunt bodies. 

Figure 26 shows the upwash calibration for several 
combinations of ports in the vertical meridian on the 
8-deg spherical nosecone at Mach 0.9. Because the 
FADS solution algorithm requires three ports on a 
vertical meridian, a total of ten different three-port 
arrangements (triples) are possible from the five ports 
shown. The variation between the corrections (8a) is 
approximately 0.5° at angles of attack less than 10°. At 
greater than 10° angle of attack, the upwash correction 
tends to be port-dependent. Port combinations that are 
distributed symmetrically about the model longitudinal 
axis (that is, triples (T, 2, 4] or [ 1,4, 6]) have corrections 
that are less than 1° to angles of attack greater than 30°. 
Triples that are clustered well below or above the 
vehicle axis of symmetry (that is, triples [2, 6, 10] or 
[1,6, 10]) have the largest corrections. 


Port 

numbers 


5a, 

deg 


2 , 6,10 



a deg 
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Initially, the calibration of the upwash and sidewash 
parameters was believed to be configuration-dependent 


Figure 26. Variation of upwash calibration with port 
selection, 8-deg cone. 
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A probable cause for the port-dependent upwash 
calibration is a variation in the vehicle-induced upwash 
from one flush-mounted port to the next. On a blunt 
body, the local upwash can significantly vary from 
location to location because the forebody is a significant 
portion of the vehicle length. With a nose diameter of 
approximately 1 in., the F/A-18 vehicle upwash cannot 
significantly vary between the ports, which would 
explain why the FADS calibration was essentially port- 
independent for this vehicle. 

An important result of the calibration dependence on 
the port locations for blunt-body shapes is that the full- 
scale vehicle must have port locations in approximately 
the same locations as the calibration model. Because the 
calibration was previously thought to be independent of 
the port locations, the flight software was thought to 
need only the location of the actual vehicle port clock 
and cone angles. This finding requires that the FADS 
ports be well-defined before the wind-tunnel calibration 
tests. Alternately, the wind-tunnel calibration model can 
be machined with additional ports (a good idea in any 
case), and the data can be used to help define the proper 
port locations on the full-scale vehicle. For the X-33 
vehicle, internal structural and space limitations also 
influenced the selection of the port locations. 

Summary 

The flush airdata sensing (FADS) system pressure 
model is derived as a synthesis of the incompressible 
potential flow around a sphere and modified Newtonian 
flow theory. The pressure model and its solution 
algorithm provide a method for extracting the entire 
airdata parameter state vector—except true airspeed— 
from flush surface pressure measurements. This paper 
uses analytical and empirical data to calibrate the FADS 
pressure model (FPM) for numerous blunt-forebody 
configurations. The sources of the forebody pressure 
data are as follows: 

• Analytical solutions for axisymmetric ellipsoids 
with fineness ratios from 1 to 20 and free-stream 
angles of attack from 0° to 70°. 

• Wind-tunnel tests of three blunt-forebody shapes— 
an 8-deg cone, a Rankine half-body, and an X-33 
forebody—with Mach numbers ranging from 0.3 to 
4.5 and angles of attack from -4° to 30°. 

• Wind-tunnel tests of X-33, X-34, and X-38 
complete-vehicle models, with Mach numbers 
ranging from 0.3 to 4.8 and angles of attack from 
-10° to 30°. 


• Flight data from full-scale tests of an F/A-18 and a 
modified F-14 aircraft. Flight test range varied from 
Mach 0.3 to Mach 2.0 and from -10° to 50° angle 
of attack. 

All of the data sets demonstrated the FADS model 
was generally applicable for blunt forebodies. The 
general structure of the calibration curves was 
consistent across the wide range of shapes, scales, Mach 
numbers, and angles of attack that were tested. No 
attempt was made to validate the model for sharp-edged 
or two-dimensional shapes. 

Specific lessons learned include the following: 

• The FPM—adapted from the solution for potential 
flow over a sphere—can be calibrated for a wide 
variety of blunt-forebody shapes. 

• The effect of lee-side flow separation is significant. 
When the surface flow has separated, the measured 
pressure cannot be accurately modeled by the FPM. 
Inclusion of pressure ports whose surface flow is 
separated will dramatically reduce the accuracy of 
the airdata estimates. Restricting the incidence- 
angle range to the ports in the attached flow region 
results in a good model of the pressure distribution. 
Using wind-tunnel data to identify the flight 
conditions at which individual ports have separated 
flow is generally conservative, because the full- 
scale flight data will usually be at higher Reynolds 
numbers (delayed separation). 

• The calibration parameters are composed of three 
terms: a shape and compressibility parameter, e; an 
upwash parameter, 5a; and a sidewash parameter, 
5(3. Using these three parameters, the FPM was 
calibrated to equivalent levels of accuracy for each 
of the data sources analyzed. When plotted as a 
function of Mach number, angle of attack, and 
angle of sideslip, the calibration parameters 
exhibited very similar characteristics. 

• All three calibration parameters were found to be 
smooth, monotonic functions that are easily 
storable in a real-time system. At high Mach 
numbers, e universally approached 0, which 
reduces the FPM to modified Newtonian flow. Also, 
for Mach numbers greater than approximately 1.5, 
the variation in e with angle of attack and angle of 
sideslip is negligible. 

• The angle-of-attack and -sideslip calibration 
parameters, 5a and 5p. are generally a function 
of angle of attack and Mach number: however, at 
high Mach numbers these values diminish to small, 
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constant-bias values. This high Mach number 
result—along with the earlier conclusion that the 
effect of angle of attack is negligible at high Mach 
numbers—greatly reduces the complexity of the 
calibration matrix, which must be tested for 
supersonic flow conditions. Furthermore, these 
results are reassuring in that they verity the FPM 
converges to modified Newtonian flow theory at 
high Mach numbers. 

• The derived calibrations are sensitive to the port 
locations. Data presented showed that the forebody 
upwash calibrations can vary by as much as 
20 percent, if ports on the test forebody do not 
match the location from which the original 
calibrations were derived. 

• Tests of complete and forebody-only X-33 
configurations revealed moderate differences in the 
calibration results at subsonic and transonic Mach 
numbers. Wind-tunnel tests of FADS calibration 
should use a complete model of the desired 
configuration. Usually, FADS calibration data can 
be obtained simultaneously with traditional force 
and moment tests. 
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